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We study fusion reactions of with ^^''Sm, ^*^W and ^-'^^U at sub-barrier energies by a coupled- 
channels framework. We focus especially on the effects of /Js deformation and low-lying vibrational 
excitations of the target nucleus. It is shown that the inclusion of Pa deformation leads to a consid- 
erable improvement of the fit to the experimental data for all of these reactions. For the ^^^Sm and 
^■^^U targets, the octupole vibration significantly affects the fusion barrier distribution. The effect 
^\ • of /3 band is negligible in all the three reactions, while the 7 band causes a non-negligible effect 

0^ ' on the barrier distribution at energies above the main fusion barrier. We compare the optimum 

0^ , values of the deformation parameters obtained by fitting the fusion data with those obtained from 

inelastic scatterings and the ground state mass calculations. We show that the channel coupling 
{^JQI of high multipolarity beyond the quadrupole coupling is dominated by the nuclear coupling and 

^ ■ hence higher order Coulomb coupling does not much influence the optimum values of Pa and Pe 



parameters. We also discuss the effect of two neutron transfer reactions on the fusion of with 

238 y 

25.7Q.Jj, 21.60. Ev, 24.10.Eq, 23.20.Js 
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\ It is now well established that nuclear intrinsic motions significantly enhance the fusion cross-section in heavy-ion 
On ' reactions at sub-barrier energies. Deformation effect is one of such prominent effects. The role of static deformation in 
On ' enhancing fusion cross section has long been recognized ||,^ and has been experimentally demonstrated [^|-^. Here 

, the enhancement occurs because there is a distribution of barrier heights which can be thought of as resulting from 

■ different orientations of the deformed target nucleus. Any distribution of barriers around a single Coulomb barrier 
I ' leads to enhancement of the fusion cross section at energies below the single barrier, because passes through the lower 

^ barriers are much more probable. Recently, high precision experimental data were obtained for the ^*^0+^*^W, "'^^'^Sm 
^ fusion reactions and it was clearly demonstrated that sub-barrier fusion reactions strongly depend upon the nuclear 
• • I hexadecapole deformation It was pointed out that the optimum values of the quadrupole and hexadecapole 

. 5^ ' deformation parameters obtained from the analyses of such high precision fusion data are consistent with those 
. obtained from the Coulomb excitation using similar radius parameter |lO| . In order to reach this conclusion, the 

■ authors of Refs. included up to the (3/^ deformation in their analyses, neglecting higher order deformations 
5^ \ such as [3q. On the other hand, the differential cross sections of inelastic alpha particle and proton scatterings 

from ^^''Sm and neutron scattering from ^^^W |l3| show important effects of [3^ deformation of the target nuclei. 
The important role of /Sg deformation of the target nucleus has been shown also in the inelastic alpha and proton 
scatterings from ^'^^U, which give optimum deformation parameters consistent to each other [p"4|-[T7|. Although each 
reaction might be sensitive to different channels, it is not obvious whether the effects of /Jg deformation on fusion 
reactions are negligible. We also notice that there still remain noticeable discrepancies between the experimental and 
theoretical barrier distributions in Ref. which included up to [3ji deformation. Although a better reproduction of 
the experimental data of the fusion excitation function as well as the fusion barrier distribution has been obtained by 
including the effects of low-lying two 2+ and one 3^ vibrations and a positive Q-value transfer channel, the role of 
higher order deformations has not yet been clarified. This motivated the present work, where we perform a detailed 
study of the effects of higher order deformation, especially the effects of (3q deformation on heavy-ion fusion reactions. 



*E-mail address: ruminOnucl. phys.tohoku.ac.jp 
^E-mail address: hagino@phys.washington.edu 
tE -mail address: takigawa@nucl.phys.tohoku.ac.jp 



1 



Besides clarifying the mechanism of heavy-ion fusion reactions, the study of the effects of higher order deformation 
is interesting to see the possibihty of heavy-ion fusion reactions as a new powerful method of nuclear spectroscopy. 
This is another motivation of the present work. We therefore compare the optimum values of the deformation 
parameters obtained from the analyses of the fusion data with those from inelastic scatterings and the ground state 
mass calculations. In order to have reliable results, one has to take various channel coupling mechanisms into account 
which might cause effects of similar order. In this connection, we discuss in this paper the effects of vibrational 
excitations of deformed targets. We also pay attention to the role of higher order Coulomb coupling. 

For these purposes, we particularly study ^^O -I- ^^'^Sm, ^^^W, ^ss^j f^gjon reactions, where high precision exper- 
imental data have been obtained We discuss the effects of channel coupling through the excitation function 
of the fusion cross section and the fusion barrier distribution [|l9|, which is defined as the second derivative of the 
product of the bombarding energy E and the fusion cross section ap with respect to E. Though the fusion barrier 
distribution has, strictly speaking, clear physical meaning only in the limit of sudden fusion, i.e. in the limit where 
the excitation energy of intrinsic excitations can be ignored, it has been shown that the concept still holds to a good 
approximation even for non-zero excitation energy [ pO[ . This method has often been used to analyse high precision 
heavy-ion fusion data and is now well known to provide a very sensitive test of various channel coupling effects. 

The paper is organized as follows. In Sec. II, we present the results of the coupled-channels analysis which takes 
only the ground state rotational band into account. The main result is that Pe deformation plays an important role in 
all the three reactions. For ^^^W, the magnitude as well as the sign of the are consistent with the results of ground 
state mass calculations and inelastic neutron scattering. However, the sign of Pq for the Sm and U targets is predicted 
to be opposite to the result of other studies. In Sec. Ill, we examine the vahdity of the calculations used in Sec. 
II, which take full order for the nuclear coupling, while only the linear order for the Coulomb coupling into account. 
By performing coupled-channels calculations keeping up to the second order in the Coulomb coupling, we show that 
the high multipolarity couplings, i.e. the Y4 and Yg couplings, are dominated by the nuclear coupling and hence 
non-linear Coulomb coupling does not almost alter the optimum values of ^4 and /Sg- In Sec. IV we present the results 
of coupled-channels analyses for ^^O + ^^^Sm and ^ssy fusion reactions which take octupole vibrations into account. 
These target nuclei have low-lying K=0~ octupole bands, which are strongly excited by the Coulomb excitation. The 
E3 transition strength from the ground state to the 3~ state is 24 W.u. for ^ssy g^j^^j^ ^ W.u. for ^^''Sm. We show 
that the octupole vibration significantly affects the fusion barrier distribution and modifies the optimum values of 
deformation parameters to fit the experimental data. Especially, it changes the sign of Pe deformation to agree with 
the analyses of inelastic a and proton scatterings and the ground state mass calculations. Notice that there is no 
experimental evidence for the low-lying octupole K=:0^ band in ^^^W, suggesting its absence in this nucleus. All 
the three target nuclei have low-lying /3 and 7 bands, whose interband E2 transition probabilities from the ground 
state 0+ to the 2+ member are: 1.0 and 4.4 W.u. in ^^'^Sm, 8.9 and 1.4 W.u. in ^^^W and 3.0 and 1.5 W.u. in 
^^^U for the 7 and f3 bands, respectively. In Sec.V, we examine the effects of /3 and 7 vibrations, and show that 
the effect of /3 band is negligible, while the 7 band affects the fusion barrier distribution at high energies. Besides 
nuclear intrinsic excitations, nucleon transfer reactions between the colliding nuclei can enhance the low-energy fusion 
cross section. In Sec. VI, we study the effect of pair neutron transfer reactions on the '^^O -l-^'^^U fusion reactions and 
discuss whether it explains the experimental fusion cross section which is systematically larger than the prediction of 
the coupled-channels calculations which ignore particle transfer reactions at low energies. We summarize the paper 
in Sec. VII. Appendix A is added to briefly explain the theoretical framework of the coupled-channels calculations we 
use. We also add Appendix B to show the structure of the higher order Coulomb coupling. 

II. EFFECT OF /Je DEFORMATION 

In this section we present the results of coupled-channels calculations which take only the ground state rotational 
band of the target nucleus into account. We treat the projectile as inert, since its excitations can be well incorporated 
with a choice of the bare potential Instead of handling the full coupled-channels equations, we introduce the 

no-Coriolis approximation throughout this paper, and ignore the change of the centrifugal potential barrier due to the 
finite multipolarity of nuclear intrinsic excitations |^,^ . This leads to considerable reduction of the dimension of 
the coupled-channels equations. We assume an axially symmetric deformation for the target nucleus and expand the 
radius up to the hexacontatetrapole deformation /Jg. We introduce the sudden tunneling approximation, and set the 
excitation energy of the ground state K=0"'' rotational band to zero. Together with the no-Coriolis approximation, 
this leads to a set of decoupled eigen-channel problems, each of which corresponds to the fusion with a fixed orientation 
of the target nucleus. Accordingly, we first solve the Schrodinger equation for a given orientation 6 for each partial 
wave J using the incoming wave boundary condition to obtain the tunneling probability Pj{E, 9). We then calculate 
the total tunneling probability Pj{E) for each J by taking average over all orientations as. 
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PAE) = IJ^ p., {E, 9) sin 9d0, (1) 

where the weight of the average has been determined by the ground state wave function of the deformed target, which 
is initially in the 0"*" state. The fusion cross section is then obtained by the standard partial wave sum. Once the fusion 
excitation function has been obtained, the fusion barrier distribution is calculated by the point difference formula of 
Ai?=2MeV in the laboratory energy, whose value was employed in Refs. [p|,[7|p0|] in analyzing the experimental data. 

We first determine the nuclear potential parameters for each target nucleus by fitting the fusion cross section larger 
than 200 mb by a potential model 0. We then calculate the fusion cross section by switching on deformations 
of different multipolarity successively. At each step, we determine the values of the deformation parameters by 
fitting of the data of fusion excitation function and readjust the potential parameters. We use the values in Ref. 

as the initial values. The results arc shown in Fig.l. The left and right columns are the excitation function 
of the fusion cross section and the fusion barrier distribution, respectively. The top, center and bottom panels are 
for ^^^Sm, i^^W and ^ssy targets, respectively. The dashed, the dotted and the sohd fines represent the results of 
coupled-channels calculations including only (32 deformation, p4 in addition and (3q as well, respectively. The dashed 
line for the i^O+^^^W fusion reactions cannot be seen clearly in the fusion excitation function, because it overlaps 
with the solid line in the semi-logarithmic plot of the present scale. However, it is clearly separated from the other two 
lines in the fusion barrier distribution. This typically shows the high sensitivity of the fusion barrier distribution to 
different channel coupling effects. The deformation parameters obtained in the analysis are shown in the figure. Those 
obtained including only up to deformation for the ^^''Sm and ^^^W targets somewhat differ from those obtained in 
[0, which are 132— 0.33, ^4—0.05 and /32=0.31, /34=— 0.03, respectively. Since we use the same radius parameter as in 
Ref. , these differences can probably be attributed to the different methods to calculate the fusion cross section in 
two works. We calculated it by numerically solving the Schrodinger equations, while Ref. introduced the parabolic 
barrier approximation, which does not work at energies far below the barrier. 

The importance of the /Jg deformation can be clearly seen in the fusion barrier distribution for all cases, and in the 
fusion excitation function as well for the ^^*Sm and ^ss^j targets. The agreement between the experimental data and 
the coupled-channels calculations concerning the fusion barrier distribution has been significantly improved above 56 
and 66 MeV for ^^^Sm and ^*^W, respectively, by including /3e deformation. The Pe deformation removes a sharp 
peak at around 82 MeV in the fusion barrier distribution, which appears in the coupled-channels calculations without 
f3e deformation, for the ^ss^j target. 

In Fig. 2, we show the dependence of the fusion excitation function and the fusion barrier distribution on the f3e 
parameter for the ^^O -t- ^^''Sm, ^^^W reactions. The three lines in each figure have been calculated by using the 
same parameter sets as in Fig.l (the solid lines), or by inverting the sign of the (3e parameter (the dotted lines) and 
by setting it to be zero (the dashed lines). This figure also shows that the effect of /3e deformation on the fusion cross 
section is not negligible. 

We compare in Table 1 the optimum values of the deformation parameters thus obtained with those obtained 
from the analyses of inelastic scatterings and the ground state mass calculations. The table also shows the radius 
parameter used in each analysis. For all three target nuclei, we observe noticeable discrepancies in the magnitudes of 
the deformation parameters among different studies (magnitude problem). This problem is, however, not so serious as 
it appears, because the discrepancies are largely due to the different choice of the radius parameter in each analysis. 
Since the strength of the channel coupling depends on the product of the deformation and radius parameters for the 
nuclear part and on f3\ x for the Coulomb part (see Eqs.(A.2-A.4)), physically important quantities are these 
products. Based on this idea. Table 1 shows the scaled deformation parameters as well, which have been calculated 
by /3\ X ro/1.06 (figures with a star) or by P\ x (rc/1.06)'^ (figures with two stars) from the original deformation 
parameters. We observe that the scaled deformation parameters from non-fusion studies are now much closer to 
the optimum deformation parameters from the fusion analysis. We cannot unfortunately rescale the deformation 
parameters for the neutron scattering from W, since the radius parameter is not given in . 

We wish to especially remark that the sign and the magnitude of (3e obtained from fusion analysis are consistent 
with those obtained from the ground state mass calculations for the -'^^^W target. Our result is consistent also with 
the neutron scattering, though it gives only the upper bound of the magnitude. On the other hand, the predicted 
sign of Pa is opposite to the results of other studies for ^^''Sm and ^ssy (sign problem). We show in Sec. IV that the 
effect of octupole vibration provides a possibility to cure this sign problem. We note that the optimum deformation 
parameters of -'^^'^Sm and ^ssy obtained from many experiments of proton and a particle scatterings agree quite well 
to each other including the sign and magnitude of (3^, though there exist a few exceptions in the case of ^^"^Sm. 

For reference, we compare in Fig. 3 the experimental data and the results of coupled-channels calculations using 
deformation parameters from the fusion data (the solid line), the ground state mass calculations ||2^ (the dotted line) 
and inelastic scatterings (the dashed line). For the latter, we choose the results of inelastic alpha [|ll|, neutron |1^ ] 
and proton ||l^ scatterings for ^^'^Sm, ^^^W and ^•^*U, respectively. In calculating the dotted and the dashed lines. 



3 



the radius of the target nucleus Rt has been adjusted such that equals that in the fusion calculations, where 

1 /3 

we chose i?T=l-06xAj, fm in order to have the same nuclear quadrupole coupling as for the solid line. Naturally, 
the parameter set obtained by the fit of the fusion data provides the best fit to the experimental data. 

Table 1: Comparison of the optimum deformation parameters and nuclear radius parameter in various analyses. 



Nuclei 


Methods 


ro (fm) 


Tc (fm) 


P2 


I3i 




^^^Sm 


^''O+^^'^Sm Fusion 


1.06 




0.322 


0.027 


0.027 




Mass Calculation l|2l 


1.16 




0.27 


0.113 


-0.005 










0.295* 


0.124* 


-0.005* 




a Scattering 


1.492 




0.225±0.005 


0.050±0.005 


-0.015±0.010 










0.317* 


0.070* 


-0.021* 




proton Scattering [f2l 






0.285 


0.051 


-0.015 




ibQ_^i8b^ Fusion 


1.06 




0.285 


-0.031 


0.027 




Mass Calculation ^ 


1.16 




0.23 


-0.107 


0.02 




neutron Scattering [p^ 






0.25* 


-0.117* 


0.022* 








0.203±0.006 


-0.057±0.006 


< 1-0.041 




WQ^AMU Fusion 


1.06 




0.289 


0.01 


0.044 




Mass Calculation jH] 


1.16 




0.215 


0.093 


-0.015 




a Scattering |0 






0.235* 


0.102* 


-0.016* 




1.2 




0.22±0.01 


0.06±0.01 


-0.012±0.01 








0.25* 


0.068* 


-0.014* 




proton Scattering [0 




1.25 


0.225±0.005 


0.045±0.005 


-0.015±0.003 










0.313** 


0.087** 


-0.040** 



III. HIGHER ORDER COULOMB COUPLING 

The results in Sec. II have been obtained by treating the Coulomb coupling in the linear order and the nuclear 
coupling in full order. Though this approximation is often used in literatures, it is worth checking the validity, 
especially in discussing the role of higher order deformations. One would guess that this approximation breaks 
down when the charge product of the projectile and target gets large. As examples, we performed coupled-channels 
calculations for the '^^S + ^^*Er and ^^0 + ^^''Sm fusion reactions by including the second order Coulomb coupling 
and by assuming only the quadrupole coupling. We found that the second order Coulomb coupling noticeably modifies 
the fusion barrier distribution. Naturally, the modification is more significant for the former reaction. An important 
issue in the context of the present paper is whether the higher order Coulomb coupling significantly changes the 
optimum values of higher order deformation parameters that reproduce the experimental data of fusion cross section. 
In this connection, we show in Appendix B the higher order terms in the Coulomb interaction up to the order of 
Pq, i.e. (34 X ^2- Although /?! would contribute in the same order, we do not show it, since it is very tedious to 
evaluate it and also its effects are negligible as we argue below. Note that other terms, like /3| and /32 x are 
higher order contributions, which are the same order of /?§ or higher, and are not shown. Eq.(B.l) indicates that 
the optimum values of /34 and pe parameters will be considerably altered by the non-linear coupling if the Coulomb 
coupling significantly contributes to the higher multipolarity, i.e. Y4 and Yg, couplings. 

In order to examine the situation, we compare in Fig. 4 the fusion barrier distribution calculated in four different 
ways. For simplicity, all the calculations have been performed by treating both the nuclear and Coulomb couplings in 
linear order and by expanding up to the Yg term. The solid line is the fusion barrier distribution obtained by keeping 
both the nuclear and Coulomb couplings as they are. The dashed line has been obtained by discarding the nuclear 
Y4 coupling term. We observe a significant change of the fusion barrier distribution. We performed additional two 
calculations, where only the Coulomb or both the nuclear and Coulomb Y4 couplings are discarded. Their results 
are almost the same as the solid and the dashed lines, respectively. These results indicate that the Y4 coupling is 
far dominated by the nuclear coupling. We checked that a similar situation holds also for the Yg coupling. We thus 
conjecture that the main effect of the Coulomb coupling resides in the Y2 coupling, and one can determine to a 
good approximation the optimum (3^ and Pe parameters through the coupled- channels calculations using the linear 
Coulomb coupling. In the following analyses, we thus treat the Coulomb coupling in the linear order. Keep, however, 
in mind that the optimum P2 value could be noticeably affected depending on whether one uses the linear or higher 
order Coulomb coupling. 
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IV. EFFECT OF OCTUPOLE VIBRATION 



We now study the effects of octupole vibration on tlie ^^O + ^^"'Sm and + ^^Sy f^jgjQ^ reactions. As already 
mentioned in the introduction, there exist low-lying K=0^ octupole bands in ^^^Sm and ^^^U, which are strongly 
excited by the Coulomb excitation through the E3 transition. We talce into account their effects on fusion by solving 
coupled-channels equations for each orientation of the deformed target. We call this procedure the 6'-scheme. We 
confirmed that the results are almost the same as those obtained by treating the rotational excitations not by the 
0-scheme, but by specifying each excited level by its spin, and by solving coupled-channels equations with a larger 
dimension which include both the K=0+ ground state and K=0^ octupole bands Similarly to the rotational 

coupling, we treat the nuclear part of the octupole coupling in full order, while the Coulomb part in the linear order. 
The amplitude of the zero point motion of the octupole vibration, which governs the strength of the channel-coupling, 
is determined from the experimental value of the reduced transition probability B(E3o T) [ p6| , p7[ from the ground 
state to the 3~ state of the K=0~ octupole vibrational band following 



^ (2) 



1 + \\J 1/^2 + |/34 + I3x?ix3 V ir/^e^ 

The optimum set of deformation parameters as well as the potential parameters are readjusted by the fitting after 
including the 3~ vibrational state. The results are shown in Fig. 5 by solid lines in comparison with the experimental 
data and the previous calculations which include only the ground state rotational band (the dashed lines). The effects 
of the octupole vibration are visible especially in the fusion barrier distribution. The resultant optimum deformation 
parameters are /32=0.314, /34=0.011 and /36=-0.016 for ^^^^Sm and /32=0.279, /34=0.0007 and /36=-0.024 for ^s^U. An 
interesting result is that the sign problem of (3q parameter has been resolved for both ^^'^Sm and ^ssy nuclei, although 
the optimum values of /34 become too small, especially in ^'^^U, compared with the other analyses. 



V. SIMULTANEOUS EFFECTS OF OCTUPOLE, p AND 7 VIBRATIONS 

We now add the effects of the (3 and 7 vibrations. We treat all the vibrational excitations by a coupled-channels 
framework by keeping their finite excitation energies and using the linear coupling approximation not only for the 
Coulomb but also for the nuclear parts. The rotational coupling is treated in the same way as in Sec. II. 

The amplitudes of the zero point motion of the (3 and 7 vibrations are determined from the experimental values of 
the reduced transition probability B(E2t) Q from the ground state to the 2+ state of the /3 band and to the band 
head of the 7 band. The formulae we use are, 

(^)(l-|v^/?2) 

Table 2 collects the experimental transition probabilities and the values of the zero point motion amplitudes for the (3 
and 7 vibrations together with those for the octupole vibration. It includes also (32 values. They have been extracted 
in Sec. II and used to determine oig , and following Eqs. |^ and ^. 

Table 2 : The zero point motion amplitude of the octupole, (3 and 7 vibrations. 



Nuclei 


i?(i?3„ T)(e26') 


B{E2p ]){e^b^) 


B{E2^ ]){e'9) 


P2 


< 




1 


^^^Sm 

238u 


0.100 |26| 
0.575 p7|l 


0.023 [ 
0.009 i 
0.0656 


28 
28 
\2l 




0.069 1 

0.150 

0.131 


28 
28 
28 


1 
1 


0.322 
0.285 
0.289 


0.103 
0.109 


0.026 
0.012 
0.0224 


0.051 
0.054 
0.034 



Since the coupling to the 7 band depends on the second Euler angle cf), we first solve the coupled-channels equations 
for a given set of (0,0) parameters. The fusion cross section for each partial wave J is then calculated by taking 
average over both and 

^^(^) = sined^y^ d<t)Pj{E,e,4>). (4) 
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The integrations are performed by Gauss quadrature. Since the numerical computation is quite heavy, we have 
not optimized the deformation parameters, but fixed them to those values used to obtain the solid line in Fig.l for 
i^O+is^W and those in Fig.5 for i^O+^^^'Sm and fusion reactions. 

The results are shown in Fig. 6. Similarly to Fig.l, the fusion excitation function and the fusion barrier distribution 
are shown on the left and right sides, respectively. The solid line for the ^^O+^^^W fusion reaction was obtained by 
including the effect of only the ground state rotational band, and is the same as the solid line in Fig.l, while the solid 
line for the -'^^O+^^'^Sm and fusion reactions was obtained by adding the octupole vibration. It is slightly 

different from the solid line in Fig.5 because of the different treatment of the nuclear coupling. 

The effect of the (3 band is very small for all three systems and invisible in the scale of Fig. 6. The dotted line was 
obtained by adding the effect of 7 band. Its effect is less important than that caused by the /34 and deformations, 
but is noticeable in the fusion barrier distribution. An interesting thing is that the 7 band does not affect the fusion 
excitation function at low energies. Consequently, its effect cannot be seen clearly in the fusion excitation function 
and concentrates in relatively high energy region in the fusion barrier distribution. This contrasts with higher order 
deformations, which affect the fusion excitation fusion at low energies as well, and hence the fusion barrier distribution 
over all energy region. 



VI. EFFECT OF PAIR NEUTRON TRANSFER 



Before we close the paper, we would like to comment on possible effects of pair neutron transfer channel on the 
fusion reactions. Refs. |29| , |30|| claim that positive Q-value pair neutron transfer channels explain the isotope effects, 
seen for example in '^^Ni + °®Ni, ^*Ni + ^^Ni and ^''Ni + ^''Ni fusion reactions by enhancing the fusion cross section 
at low energies in ^*Ni + ^^Ni collision. Similarly, by studying ^*Si + ^*Zn scattering at sub-barrier energies, Ref. 
|^l| claims that the coupling of the positive Q-value two neutron transfer channel significantly enhances the fusion 
cross section. 

Among the three reactions which we discuss in this paper, only the ^^O -I- ^38^ j^^^g g. neutron transfer channel 
whose Q-value is positive, the Q-value for the two neutron pick-up reaction from ^^O -|- ^asy isq _|_ 236y being 
0.826 MeV in the ground state channel. This transfer channel might resolve the discrepancy between the experimental 
data and the coupled-channels calculations in the fusion excitation function at low energies for the ^^O -I- ^38^ fusion 
reaction (see Fig. 6). In order to see this possibility, we study here the effects of this transfer channel following the 
prescription in Ref. pof , where the transfer reaction is treated in the same way as a vibrational excitation in the 
coupled-channels formalism. The form factor of the transfer reaction is assumed to be 

Ftrans[R,0) = -at — , (5) 

an 

where at is the strength parameter of the transfer reaction, Vn{R,0) is the deformed ion-ion potential. This form 
factor is slightly simplified from that in Refs. p^ , ^ by ignoring a small correction term. We determine the strength 
parameter by fitting the excitation function of the fusion cross sections. The optimum set of deformation parameters 
are readjusted by the fitting to the experimental data after including pair neutron transfer. 

The results are shown in Fig. 7. The solid line includes the effect of the transfer reaction, while the dashed line 
takes only the rotational excitation into account. We left out the effects of vibrational coupling in these calculations. 
We see that the transfer channel significantly enhances the fusion cross section at low energies. The optimum de- 
formation parameters in this analysis are /32=0.299, P4 — 0.002 and /36=0.034. Unfortunately, the dimension of the 
coupled-channels calculations becomes too large to simultaneously take both effects of transfer reaction and vibra- 
tional excitations into account. In the summary section, we also mention the possible importance of single nucleon 
transfer reactions. 



VII. SUMMARY 



We studied the effects of Pe deformation on heavy ion fusion reactions at energies near and below the Coulomb 
barrier by analyzing the excitation function of the fusion cross section and fusion barrier distribution for ^^0+^^^Sm, 
186^ and 238^ j-eactions. Coupled-channels equations have been solved by taking the rotational excitation, and the 
octupole, /3 and 7 vibrations of the target nuclei into account stepwise. The calculations which took only the ground 
state rotational band into account showed that Pe deformation is important for all three reactions. The optimum 
value of Pe well agrees with that obtained from the ground state mass calculations and inelastic neutron scattering 
for ^^^W. On the other hand, the sign of /3e is inconsistent with that from the analyses of inelastic alpha and proton 
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scatterings and the ground state mass calculations for -'^^'^Sm and ^ssy targets. We examined the validity of the 
linear approximation we took for the Coulomb coupling and gave a reasoning to conjecture that it is good enough to 
determine the optimum values of higher order deformation, i.e. f]^ and /Jg, parameters. 

We have then shown that the coupling to the low-lying octupole vibration significantly affects the fusion barrier 
distribution in the ^®0+^^'*Sm, ^38^ reactions. Interestingly, it changed the sign of the optimum /Jg to agree with that 
suggested from non-fusion analyses. The (3 and 7 vibrations are then also taken into account. We found that the P 
vibration introduces only negligible effect, while the 7 vibration changes the fusion barrier distribution by a noticeable 
amount, though the change is less than that due to higher order, i.e. P4 and Pe, deformations. An interesting feature 
is that the 7 band does not affect the fusion excitation function at low energies, and hence its effect mostly appears in 
the fusion barrier distribution in relatively high energy region. This contrasts with higher order deformations, which 
infiuence the fusion cross section at low energies as well, and hence the fusion barrier distribution over all energy 
region. We left fitting to optimize the deformation parameters to a future work because of the computational 
heaviness. 

A problem with the fusion reaction is that the coupled-channels calculations which include only rotational 

and vibrational excitations cannot reproduce large experimental fusion cross section at low energies. We showed in 
sect. VII that two neutron transfer reaction enhances the fusion cross section at low energies. One will, however, need 
to study the effects of single nucleon transfer reactions as well in order to reach a comprehensive understanding of this 
reaction and to draw a conclusive conclusion. In this connection, it is interesting to notice that larger experimental 
cross sections for one nucleon transfer reactions than those for two nucleon transfer reactions at low energies have 
been reported for several systems 1 34 - 3^] . 

Finally, we wish to make some comments on the limitations of our theoretical framework. We assumed a simple 
Coulomb interaction given by Eq.(A4), which has a few shortcomings. The first is that the bare Coulomb interaction 
is identified with the Coulomb interaction between two point charges instead of the Coulomb potential for a uniformly 
charged extended object, which is often used for heavy-ion collisions. The second is that the same Coulomb coupling 
form factor, which is valid only in the region where there is no overlap between the projectile and target nuclei, is used 
over all separation distance. Furthermore, the Coulomb and the nuclear deformation parameters are assumed to be 
the same. The first two problems are related to each other, and might cause non-negligible effects, since the turning 
point at the inner side of the potential barrier is located sometimes inside the standard choice of the so called Coulomb 
radius, Rc = 1.2(^p ^ -I- A!j^^) fm, or even the sum of the radii of the projectile and target nuclei. The analyses which 
allow different values for the Coulomb and the nuclear deformation parameters carry one of the important advantages 
of heavy-ion fusion reactions compared with the other analyses, say neutron scattering. Such analyses will explore the 
difference between the charge and matter distributions, and will be very interesting also in connection with the study 
of the structure of exotic unstable nuclei, which is one of the current interests of nuclear physics. We will discuss in 
detail the eff'ects of the improvements of theoretical analyses in these three respects in a separate paper [|o| . Here we 
simply wish to mention that the first effect, i.e. the difference of the Coulomb interaction between two point charges 
from the Coulomb potential for a uniform charge distribution, can be mimiced to a large extent by renormalizing the 
bare potential. 
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APPENDIX A: BRIEF DESCRIPTION OF THE COUPLED-CHANNELS FORMALISM 

In this appendix we briefiy explain the coupled-channels formalism which we used. We present here the case, where 
all the /3, 7 and octupole vibrations are taken into account. The total Hamiltonian reads 

H = T + H,nt{0 + ViR,0, (Al) 

where T is the kinetic energy of the relative motion between the projectile and target, -ffmt(C) the Hamiltonian of the 
intrinsic motions of the colliding nuclei, whose coordinates are denoted by ^, and V{R, ^) the interaction Hamiltonian 
which depends on the coordinates of both the relative motion, R, and nuclear intrinsic motions. 

We use the geometrical collective model for nuclear intrinsic motions. The variables ^ are then the static as well as 
dynamic deformation parameters specifying the radius of the target nucleus as. 
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R{e, (j), a) = Rt 



1 + PxyXoiO) + 020^20 {0) + ^22 (0, 0) + F2-2 {0, 0)] + 030^30 {0) 



(A2) 



In writing Eq.( A2) we chose the rotating coordinate frame, where the z-axis is taken to be paraUel to the coordinate of 
the relative motion R |2^ ] . 9 and (f> are Euler angles which define the orientation of the principal axes of the deformed 
target in this frame. P\, A being 2, 4 and 6, are the static deformation parameters, a means a2Q, and 030 which 
are the dynamical deformation parameters describing the /3, 7 and octupole K=0~ vibrations, respectively. 
The interaction Hamiltonian V{R,£,) consists of the nuclear and Coulomb parts. We assume the former to be 



Vn{R,0,(I>, a) 



1 + exp[(i? ~ Rp- R{e, </), a))/ao] ' 



(A3) 



It contains both the bare potential and the coupling Hamiltonian. When we treat the channel-coupling in the 
perturbation theory, say, of first or second order, we expand Vn{R, d, 0, a) with the relevant deformation parameters. 
The actual procedure of the full order coupled-channels calculations is explained in Refs. |37|j3^ ]. 

We assume a simple Coulomb interaction by ignoring the change of the analytic expressions of the bare Coulomb 
interaction and the Coulomb coupling form factor depending on the relative magnitude between the distance R and 
either the sum of the charge radii of the projectile and target or the absolute value of their difference. The formula 
we take reads up to the leading order of the dynamical variables as 



Vc{R.e,4>,a) = 



ZpZre^ 
R 



E 



2A 



"iZ pZj'C^ R^ 
~R>^ 

Yx-2 



+Sx,2a'x2 (^2 (9, 0) + Yx-2 [0, 0)) 



^ZpZre'^ ,„.Rt 
^ asoYsoiO)^. 



(A4) 



We assume the same charge radius and deformation parameters as those for the nuclear part for the target nucleus. 
The Hamiltonian for the intrinsic motions consists of four parts. 



Hint {£.) — Hrot + Hp + Hy + Ho 



(A5) 



They describe the rotational and /3, 7 and octupole vibrational excitations. Their explicit forms and the corresponding 
eigenfunctions and eigenvalues can be found in Ref. 

We introduce two basic approximations. The one is the no-Coriolis approximation and the other the sudden 
tunneling approximation, i.e. degenerate spectrum approximation, for the rotational motion. The latter corresponds 
to setting Hrot to be zero. In these approximations, coupled-channels equations are solved for each given set of 
(J, 0,0), J being the initial angular momentum of the relative motion, by expanding the wave function as 



R 



*n(^„), 



(A6) 



where ^t, represent the coordinates of the /?, 7 and octupole vibrations, and n is the abbreviation of a set of corre- 
sponding quantum numbers {np, n-y, Uo) [ p9[ . Note that the angular part of the wave function for the relative motion 
is simply a constant in the rotating frame approximation. We consider only or 1 for all the vibrational quantum 
numbers. The coupled-channels equations read 



2fi dR2 



J{J + l)h^ 



+ em + V{R,e,a^O) ~ E 



Vmn{R,9,^)xi''^{R) 



(A7) 



where is the eigenvalue of the vibration excitations corresponding to the eigenstate ^m{£.v)- We represented the 
total interaction by separating it into the diagonal V{R,9,a — 0) and the explicit coupling Vmn terms with respect 
to the vibrational excitations. The latter have been evaluated using the wave functions for vibrational motions given 
in Ref. ||. 

We solve the coupled-channels equations by imposing the incoming wave boundary condition at the position of the 
s-wave potential minimum, and determine the fusion probability by evaluating the incoming flux in each channel at 
that position. Once the fusion probability is obtained in this way for a given set of (J, 0,(j}), the total fusion probability 
for that partial wave is calculated by taking average over the orientation (0,0) as given by Eq.(|). The fusion cross 
section is then calculated by the usual partial wave sum. 
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APPENDIX B: HIGHER ORDER COULOMB COUPLING 



Here we present the explicit form of the Coulomb coupling up to the Yg term when the second order coupling terms 
are included. Only the major terms are explicitly shown for the second order coupling. 

^ ( n a ^ 20 /45 X 13\ ^ ^ 2^T,^ 
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Figure Captions 

Fig.l 

Comparison of theoretical (a) excitation functions of the fusion cross section and (b) fusion barrier distributions with 
experimental data [[zljl^]. Only ground state rotational excitations are taken into account in the theoretical calcula- 
tions. The dashed lines represent the optimum fits when only quadrupole deformation is included. The dotted curves 
include hexadecapole deformation, while the solid lines show the final fits including hexacontatetrapole deformation. 

Fig.2 

Sensitivity of (a) the excitation function of the fusion cross section, and (b) the fusion barrier distribution to the 
hexacontatetrapole deformation. The solid lines were obtained by using the optimum /Jg parameter. The dotted lines 
were obtained by inverting the sign of /Jg, while the dashed lines represent the results when the hexacontatetrapole 
deformation is set equal to zero. 

Fig.3 

Comparison of theoretical (a) excitation functions and (b) fusion barrier distributions calculated by using deformation 
paramaters from different analyses. The data are taken from Refs. |0Jl^. The results using the optimum deformation 
parameters from the analysis of the fusion data are shown by the solid lines. The dashed lines use the deformation 
parameters obtained from inelastic scattering |0,0,n^ , while the dotted curves those from the nuclear ground state 



mass calculations 24 1. All the calculations include hexacontatetrapole deformation. 



Fig.4 

Study of the relative importance between the nuclear and the Coulomb I4 couplings. The solid line has been calcu- 
lated by including both nuclear and Coulomb Y4 couplings, while the dashed line by ignoring the nuclear Y4 coupling. 
The results, where only the Coulomb Y4 coupling and both the nuclear and Coulomb Y4 cuplings have been discarded 
are almost the same as the solid and the dashed lines, respectively. 

Fig.5 

Effects of octupole vibration. The dashed lines represent the optimum fits when only the ground state rotation band 
is taken into account, while the solid lines show the optimum fits when the coupling to the K=0~ octupole vibration 
is added. 



Fig 6 

The same as Fig.l, but when octupole, /3 and 7 vibrations are taken into account. The solid line for the ^^^W target 
has been obtained by including only the ground state rotational band, while that for the ^^''Sm and ^ssy targets by 
adding the octupole vibration. The change due to the beta band is invisible. The dotted line has been obtained by 
adding the 7 band. 

Fig 7 

Effect of pair neutron transfer reaction on the ^^0-|-^'^^U fusion reactions. The dashed line takes only the ground 
state rotational excitation into account, while the solid line the pair neutron transfer reaction in addition. 
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